function y = XtrackMDLpara(model)

    % robust edge neighbor threshold 
    d = 10;
    
    model_para = zeros(3,size(model,3));

    for i=1:size(model,3),
        
        Ax = model(1,1,i);
        Ay = model(2,1,i);
        Bx = model(1,2,i);
        By = model(2,2,i);
        
        dx = Bx-Ax;
        if (dx~=0),
            a = (By-Ay)/dx;
            b = (Bx*Ay-Ax*By)/dx;
            eps = d*sqrt(a^2+1);
        else
            a = inf;
            b = 0;
            eps = d;
        end
        
        model_para(:,i) = [a b eps]';
    end
       
y = model_para;

end